Oil displacement through a porous medium with a temperature gradient 
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We investigate the effect of a temperature gradient on oil recovery in a two-dimensional pore- 
network model. The oil viscosity depends on temperature as, /i = exp(B/T), where B is a physico- 
chemical parameter depending on the type of oil, and T is the temperature. A temperature gradient 
is applied across the medium in the flow direction. Initially, the porous medium is saturated with 
oil and, then, another fluid is injected. We have considered two cases representing different injection 
strategies. In the first case, the invading fluid viscosity is constant (finite viscosity ratio) while in 
the second one, the invading fluid is inviscid (infinite viscosity ratio). Our results show that, for 
the case of finite viscosity ratio, recovery increases with AT independently on strength or sign of 
the gradient. For an infinite viscosity ratio, a positive temperature gradient is necessary to enhance 
recovery. Moreover, we show that, for AT > 0, the percentage of oil recovery generally decreases 
(increases) with B for a finite (infinite) viscosity ratio. Finally, we also extend our results for infinite 
viscosity ratio to a three-dimensional porous media geometry. 

PACS numbers: 



I. INTRODUCTION 

Thermal recovery processes have been used by 
petroleum companies as a strategic method to improve 
oil production from reservoirs. These processes consist, 
basically, in decreasing the oil viscosity (increasing the 
pressure) by increasing the temperature of the reservoir 
using a heat source In practice, this can be done by 
injecting a hot fluid (steam or water) into the reservoir. 
A method often called Steam or Hot- Water Injection has 
been mostly used by companies exploiting heavy oil reser- 
voirs nij. 

The temperature dependence of the oil viscosity is gov- 
erned by physico-chemical parameters which can allow 
a reduction of oil viscosity by several orders of magni- 
tude with only a modest increase of temperature In 
general, oil properties in the operational conditions of a 
reservoir field are very difficult to predict and can be very 
different for each type of oil or oil mixture. Therefore, 
their direct measurement is highly desirable for a better 
understanding of this phenomenon, and to improve the 
efficiency of the recovery. 

Several approaches to model and simulate oil recovery 
have been utilized in the past. Some of them make use 
of the macroscopic description of conservation laws in a 
porous medium, simulating a whole reservoir, including 
injectors and producers wells 0, Q. Other authors use 
conservation laws under a more microscopic approach 0- 
Il8j , where a portion of the porous medium can be repre- 
sented by tubes connected to one another. The fluid flow 
in each tube is easily computed by the Hagen-Pouseille 
equation. For instance, Lu et al. |19l - [2l| have used such 
type of modeling approach to simulate oil burning by 
air injection. This represents another thermal recovery 
method, where the heat source is the burning oil itself. 
In their work, they have studied the penetration of the 
burning front into a solid oil phase. In the present work, 



we study the displacement in a porous medium of an oil 
with temperature-dependent viscosity being pushed at 
microscopic scale by other fluid. We adopt a simple two- 
dimensional network model [13, HH previously developed 
to simulate two-phase flow with arbitrary viscosity ratio. 
Despite its simplicity, the model is capable to reproduce 
a large variety of experimental results [22j . In order to 
adapt this model to our purpose, we implement a tem- 
perature gradient in the injection direction and assume 
that the oil viscosity has an exponential dependence of 
on the inverse of temperature We studied two dif- 
ferent cases according to the viscosity ratio. In the first 
one, the invading fluid viscosity is constant (finite vis- 
cosity ratio) and in the second one, the invading fluid is 
inviscid (infinite viscosity ratio) . The aim of this work is 
then to investigate the influence of a temperature gradi- 
ent on the efficiency of oil recovery under these different 
conditions. 



II. MODEL FORMULATION 

The disordered porous medium is represented by links 
and nodes. The links represent pieces of rock of equal 
length i and cross-section area a, and to each one we 
assign a permeability k, which is chosen randomly ac- 
cording to a uniform distribution in the interval [10~ 5 , 1]. 
This randomness in the permeability represents the dis- 
order of the porous medium. The nodes where four links 
meet are assumed to have no volume. The links are 
placed on a square lattice tilted by 45 degrees which 
assures that all links are geometrically equivalent with 
regard to the average flow, i.e., the links are neither par- 
allel nor perpendicular to the flow direction. Initially, the 
porous medium is fully saturated with oil and periodic 
boundary conditions are applied at the top and bottom 
of the system. The penetration process starts with an in- 
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vading fluid being injected at constant now rate through 
the left boundary of the system. 

The volumetric flow rate in a link connecting neighbor 
nodes i and j is given by the Darcy's law, 



where pi is the pressure at node i and fey is the perme- 
ability of the link. The effective viscosity /i e f for a given 
volume link is calculated according to a linear mixing 
rule, namely, /x e / = S ^ + Sj/ij, where S and Si are 
the saturations of the oil and the invading fluid, respec- 
tively, and fj, and fii are their corresponding viscosities. 
For each species, the saturation is calculated at a link as 
the volume fraction of the corresponding phase. We as- 
sume that a link is immediately accessible when touched 
by the invading fluid at one end, neglecting wetting or 
drying effects, the pinning of the interface due to impu- 
rities, and finite contact angles or surface tension at the 
pore level. The capillary forces in our model are implicit 
in the permeability of the link and depend thus on the 
scale of this link. Therefore we cannot calculate an ex- 
plicit flow rate in physical units. 

Mass conservation at each node of the lattice leads to 
the following set of coupled linear algebraic equations: 

5^?«=0, for i = l,2,...,N, (2) 
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where N — L 2 is the total number of nodes, L is the 
linear size of the lattice in the x-direction, and the sum- 
mation j runs over the nearest neighbor nodes of node i. 
These equations are solved to obtain the node pressures 
at each time step. In order to simulate the dynamics of 
viscous invasion, we neglect the effects of fingers in each 
link containing both phases and consider an abrupt sat- 
uration profile along its axial direction. Thus we define 
an interface inside the link separating the oil and the in- 
vading fluid that is an approximation and should not be 
confused with a meniscus in a pore. Then, we allow those 
interfaces to displace by a length, Axy = qijAt m i n /a, 
where At m i„ is the minimum time, among all links con- 
taining both phases, necessary for the invading fluid to 
reach the end of a link. When an interface reaches the end 
of a link, reaching a point-like node, it is instantaneously 
transferred to those neighbour links whose pressure dif- 
ferences allow for oil displacement. To avoid multiple 
interfaces in a single link the following rule is adopted. 
When a third interface appears in a link, these three in- 
terfaces are reduced to a single one by merging bubbles 
of the same phase, so that the phase volume in each link 
is conserved. The unphysical jumps on the pressure re- 
sulting from this reorganization scheme represent only 
negligible perturbations due to the small size of the cor- 
responding bubbles [22]. This procedure is executed at 
each time step until breakthrough happens, i.e., the in- 
vading fluid just reaches the other end of the system . 
Here we assume that the viscosity of the fluid phases 




FIG. 1: Snapshots of the invading fluid near breakthrough for 
the finite viscosity ratio case, using distinct types of oil and 
temperature differences. All patterns have been generated 
with the same distribution of permeability for the random 
links and L — 80. From top to bottom, B changes from 7 to 
3, while from left to right, AT assumes the values -4, and 
4. 



typically obey an exponential dependency on the inverse 
of temperature, 

/i D = exp(B/T), (3) 

where the controlling physico-chemical parameter B, in 
units of temperature, is considered to be constant. One 
calls a fluid "heavy oil" if B is high and "light oil" if 
B is small. In order to impose a temperature gradient 
across the medium, different temperatures are assigned 
to the left (inlet) right (outlet) boundaries of the lat- 
tice and we consider, then, a linear temperature variation 
from inlet to outlet. This imposed gradient is constant 
in time, which is analogous to assume that the thermal 
conductivities of both fluids are negligible compared to 
that of the rock. The temperature difference is defined 
as AT = Ti n i e t — T out iet, and we carry out simulations 
many different dimensionless values of Ti n \ e t and T OMt ; et 
{Tirdet, Toutiet = 1, 2, 3, 4 and 5) such that -4 < AT < 4. 
Negative values of AT represent a cold injection which 
is not of technological interest, but can help to under- 
stand the nature beyond the fluid displacement under a 
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FIG. 2: Near-breakthrough saturation of the invading fluid 
for the finite viscosity ratio case. The curves correspond to 
three values of the temperature difference, AT = —4, 0, and 
4, B = 5 and L = 80. In the inset, we also show the near- 
breakthrough saturation behavior, but for a fixed value of 
AT = 4, and B = 3, 5 and 7. 



temperature gradient. In this model, both link length, t, 
and reference temperature, T r , are adjusting parameters 
which one can use to apply that model to a practical ap- 
plication. For example, if one considers T r = 15°C and 
t = 1cm (size of the fraction of the rock which each link 
represent in the model), the temperature gradient varies 
from to 75 °C/m as AT changes from to 4. In the 
next section, we show results for the two different cases of 
the viscosity ratio. In both, \x Q is given by Eq. (3) while 
\xi is constant, i.e., independent of the temperature. This 
can be justified since, for example, water, which is a typ- 
ical invading fluid, has a viscosity that hardly changes, 
compared to many types of oils, for common operational 
temperature intervals. 



III. RESULTS 

In this section, we study the patterns of the invaded 
region and the percentage of recovered oil when an exter- 
nal fluid with unity viscosity is injected into the medium. 
According to our model, the viscosity ratio, \i j has 
a finite value which changes on x-direction, if AT is dif- 
ferent from zero. In Fig. [T] we show six different patterns 
of the invading fluid with the same distribution of ran- 
dom links for different values of AT and B. From top 
to bottom we change B from 7 to 3, while from left to 
right AT assumes the corresponding values of —4, and 
4. We clearly observe that either decreasing B or for pos- 
itive values of AT, the invading patterns become more 
compact. 

In the isothermal case, AT = 0, the oil viscosity is a 
constant depending only on B. For a heavy oil (B = 
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FIG. 3: Percentage of recovered oil versus AT for the finite 
viscosity ratio case using three different values of B and L = 
80. 



7), we can see finger patterns appearing in the invaded 
region. When AT > 0, the viscosity ratio is small on the 
left side of the lattice and becomes larger as the invading 
fluid penetrates in x-direction. First, the front advances 
compactly but at some time, a finger appears and grows 
faster until it reaches the other end of the lattice. In 
the opposite case, when AT < 0, fingers appear initially 
and then become broader. For B — 3, the changes in the 
viscosity ratio due to the temperature difference suppress 
the appearance of fingers. 

For each set of parameters, we perform simulations 
with 50 realizations of the disordered porous media, to 
obtain average values of both the invading fluid satura- 
tion and the percentage of recovered oil at the break- 
through. In Fig. [21 we show the near-breakthrough satu- 
ration profile of the invading fluid for different values of 
AT and B. We observe that for a positive AT the sat- 
uration profile has a plateau before it starts decreasing 
in the x-direction. However, for a negative AT, finger 
patterns in the beginning of the lattice produce a dip in 
the saturation profile close to the inlet. The percentage 
of recovered oil versus AT is shown in Fig. [3] for three 
different values of B. We see that, despite heavy oil re- 
covery is smaller, all types of oil tend to the recovery 
performance for high values of AT. 

Now, we study the recovery of an oil that is much 
heavier than the invading fluid, i.e., for the case of an 
extremely large viscosity ratio, /x //i/ — > oc. This ideal- 
ized condition is implemented here by considering that 
the pressure in the invaded fluid immediately adjusts to 
the injection pressure, i.e., the invading fluid is inviscid. 
This simplification allows us to simulate bigger lattice 
sizes, since we need to solve Eqs. (2) only for non-invaded 
sites. 

In this case, the observed patterns of the interface 
are always viscous fingering-like, as shown in Fig. 21 for 
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FIG. 4: Snapshots of the invading fluid near breakthrough 
for the infinite viscosity ratio case. All patterns have been 
generated with the same distribution of permeability for the 
random links, B — 5 and L = 256. From left to right, AT = 
-4, and 4. 
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FIG. 5: Near-breakthrough saturation of the invading fluid 
for the infinite viscosity ratio case. The curves correspond 
to three values of the temperature difference, AT = —4, 0, 
and 4, B — 5 and L — 256. In the inset, we also show the 
near-breakthrough saturation behavior, but for a fixed value 
of AT = 4, and 5 = 1,3 and 5. 



AT = -4, and 4, with B = 5 and L = 256. For 
AT = 0, these patterns show fingers which agree with 
well know two- pha se displacement patterns with infinite 
viscosity ratio [24 l25j . We can also see that different 
patterns occur for different values of AT. The reason for 
this behavior is that, despite the inviscid characteristic 
of the defending fluid, the oil viscosity has a finite value 
given by Eq. (3) which falls as the temperature is raised. 
Then the number of longer fingers for B = 5 increases 
with AT as shown Fig. [U This observation is also valid 
for B = 1 and 3 (not shown), and therefore represents a 
standard behavior in the case of infinite viscosity ratio. 




FIG. 6: Percentage of recovered oil versus AT calculated in 
the infinite viscosity ratio case, for different values of L and 
B. In the main plot, after proper rescaling, we show that 
these curves collapse on the top of each other and can be all 
fitted to a hyperbolic tangent function. In the inset, we show 
the same results before rescaling. 



For each set of parameters, we averaged over 100 real- 
izations to obtain the invading fluid saturation and the 
percentage of recovered oil. The near-breakthrough sat- 
uration profile of the invading fluid is shown in Fig. [5] for 
different values of the relevant parameters. For a neg- 
ative value of AT, the saturation profile always decays 
in x-direction. However, for a positive AT, we can iden- 
tify a region in the center of the medium with approxi- 
mately constant saturation. The extension of this region 
increases with B because heavier oils create slower and 
wider fingers, which tend to make the invasion patterns 
more compact. 

In the inset of Fig. [SJ the percentage of recovered oil 
R{%) versus AT is shown for different values of L and B. 
The main plot shows that, after rescaling, these curves 
collapse on the top of each other and can be closely de- 
scribed by a hyperbolic tangent function in the form, 



R{%) 



a + 65 tanh(AT) 



(4) 



with parameters a — 51.11 and b — 1.74 obtained 
through the best nonlinear fit to the data. The expo- 
nent a is found to be about 0.2 and can be computed as 
a = d—df, where d = 2 is the Euclidian dimension of the 
lattice and df is the fractal dimension, which is found to 
be 1.8 for our system. The results show that above a cer- 
tain value of temperature difference, AT « 2, no relevant 
change in R(%) is observable anymore for all types of oil. 
This means that a too strong temperature gradient can 
be an unnecessary cost to the recovery process. 

Finally, we also performed simulations with our model 
in the infinite viscosity ratio regime using a more real- 
istic three-dimensional porous medium substrate. More 
precisely, we obtained results for a cubic lattice with size 
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FIG. 7: Near-breakthrough saturation profile of the invad- 
ing fluid penetrating in a cubic lattice for the case of infinite 
viscosity ratio. The curves are the results of simulations for 
AT = —4, 0, and 4, with parameter B — 5, and system size 
L = 20. In the inset, we also show the near-breakthrough 
saturation behavior, but for a fixed value of AT = 4, and 
B = 3, 5 and 7. 



L = 20 and averaged over 100 realizations for several 
values of B and AT. As shown in Figs. [7] and the sat- 
uration profiles and the recovery performance of the sys- 
tem, respectively, remain basically the same as compared 
to the results found for two-dimensional porous media 
models (for comparison, see Figures [5] and [6]) . This qual- 
itative similarity reinforces the validity of our approach 
as a way to increase the efficiency of the recovery process 
by means of a temperature gradient. 

IV. CONCLUSIONS 

In summary, the main purpose here was to investigate 
how the front between two immiscible fluids propagates 
in a model porous medium as a function of the ratio 
of their viscosity and under the influence of an imposed 
global temperature gradient. An obvious technological 
application of our study would be to enhance the re- 
covery efficiency of oil being pushed by hot-water in a 
petroleum reservoir. In our simulations, this has been 
accomplished by explicitly coupling the oil viscosity with 
the inverse of temperature locally in terms of a simple 
exponential dependence. We thus proceeded with the 
displacement of several types of oil through many real- 
izations of disordered porous medium, and subjected to 
a range of distinct temperature differences. 

Two different regimes of viscosity ratio have been stud- 
ied. In the first, oil viscosity changes with temperature 
and the viscosity ratio is always finite, but can vary over 
several orders of magnitude. In the second regime, oil is 
assumed to be "heavy" , in the sense that it is a extremely 
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FIG. 8: Dependence of the percentage of recovered oil on the 
temperature difference AT or different values of B on a cubic 
lattice with L — 20 for infinite viscosity ratio. 



viscous when compared do the invading fluid, even if a 
maximum temperature difference is applied, hence the 
viscosity ratio can be considered as infinite. We find that 
the best conditions for recovery are significantly depen- 
dent on the adopted regime. In the finite viscosity ratio 
case, an oil with a viscosity that is only weakly dependent 
on the temperature is better recovered, independently on 
the strength or sign of the gradient. Also in this case, dif- 
ferent invasion patterns can be observed as the viscosity 
ratio changes, namely, we find fronts that are compact, 
unstable and sometimes a mixing of both. 

Differently, in the infinite viscosity ratio case, oil re- 
covery increases with the exponential parameter B if the 
temperature difference is positive. Moreover, recovery is 
found to follow a hyperbolic tangent behaviour on the 
temperature difference and the best recovery is obtained 
for positive temperature gradients (Hot- Water Injection). 
It would be interesting to verify our predictions with ex- 
perimental results. Since field experiments in this area 
are usually difficult and expensive, it would be advisable 
to run laboratory-size experiments in which interface pat- 
terns are dynamically registered and the recovered vol- 
ume is systematically measured while a viscous fluid is 
displaced by a less viscous one in a porous medium under 
the influence of temperature gradients. 
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